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Abstract 



We developed a new method based on functional integration to treat the 
dynamics of polarons in one-dimensional systems. We treat the acoustical 
and the optical case in an unified manner, showing their differences and sim- 
ilarities. The mobility and the diffusion coefficients are calculated in the 
Markovian approximation in the strong coupling limit. 
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1 Introduction 



The main task of this paper is to develop a new formahsm to treat the dynamics of 
acoustical or optical polarons in the strong coupling limit. Although we are applying 
the formalism to a specific problem we think that it can also be applied to a large 
class of phenomena, specially problems involving quantization of zero frequency 
modes in theories which have solitons as solutions of their semiclassical equations of 
motion. 

When we have a particle (electron) interacting with a given background (phonons 
in this case) and wish to study its effective dynamics, it is now well-known that 
we must trace over the phonon coordinates and study the time evolution of the 
reduced density operator of the electronic system. However, in the specific case of 
electron-phonon coupling, the "effective propagator" for this operator is extremely 
cumbersome, preventing us to get any simple result out of this standard analysis. 
Therefore, one should search for an extra step before blindly tracing the phonon 
coordinates out of the problem. 

We start from a very intuitive point, treating the electron-phonon Hamiltonian 
as a "semiclassical" Hamiltonian [Q. This "semiclassical" picture provides us with 
solutions which are solitons, that is, solutions which do not change their shape with 
time. These solitons will be the basic entities for the future solution of the problem. 
We show that the best basis in which one can expand the field operators of the 
quantum Hamiltonian is obtained from the problem of an electron trapped in a 
self consistent potential well. We will call it the "adiabatic basis" since the strong 
coupling limit is the adiabatic limit (see section 2). 

Once we have obtained the Hamiltonian in the adiabatic form we can eliminate 
the electronic part perturbatively, that is, we trace over the electron coordinates. 
This treatment gives rise to an effective Hamiltonian for the phonon system which 
has renormalized phonons and a zero frequency mode, namely, the polaron. An 
important feature of this Hamiltonian is that it can be straightforwardly generalized 
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for a non-interacting many electron system . 

Using the well-known "collective coordinate formalism" we transform the effec- 
tive Hamiltonian into a Hamiltonian of a particle, the polaron, coupled to a new set 
of phonons. It can be shown that for a small polaron momentum the problem can 
be put in a very simple form. 

Actually, this is a more systematic way to apply the ideas used by Schuttler 
and Holstein to a quantum dissipation problem as the necessary step before the 
tracing procedure. Here, we shall repeat part of their arguments for the sake of 
completeness. 

Finally, using the functional integral formalism, we can show that the polaron 
behaves as a brownian particle due to the scattering of phonons. So, we have devel- 
oped a method which allows us to calculate the physical quantities of interest, such 
as the damping parameter (mobility) and the diffusion coefficient without appealling 
to kinetic theory. 

In section II we will present the model and exhibit the adiabatic basis for the 
strong coupling limit while in section III we obtain the effective Hamiltonian for the 
polaron coupled to the renormalized phonons. In section IV we use the functional 
integral formalism in order to show how this problem can be treated as a brownian 
motion problem and in section V we use the previous results to calculate the physical 
quantities of interest. Section VI contains our conclusions. 

2 The Polaron Model and the Adiabatic Expan- 
sion 

Since in this paper we will treat the problem of an electron coupled to acoustical 
or optical lattice vibrations (phonons) in an unified manner, we decided to develop 
these two problems in parallel, in order to show their differences and similarities. 
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2.1 Optical Case 

The optical polaron model is based on the Frohlich Hamiltonian [Q] for electrons 
coupled to longitudinal optical phonons. This Hamiltonian can be written in the 
second quantized form, in one dimension, as 



_L n nn' 1 



here vr and fj are the momentum and position operator for the phonon field and ip'^ 
and "ip the creation and destruction operators, respectively, for the electron field. 
They obey the following commutation rules 

[i7(x, t), vr(x', if:)] = ih5{x — x') 

{V^(a;,t),?/'^(x',t)} = 5{x - x') 

where [ , ] denotes commutation and {, } anticommutation. All the other commu- 
tation (or anticommutation) relations are zero. 

In (2.1.1), ujo is the frequency of the phonons, v = M/a is the lattice density, 
M the ion mass, a the lattice parameter, m is the effective mass for the electrons in 
the conduction band and D the coupling constant. 

In order to analyze the physical content of (2.1.1) we will treat the operators as 
ordinary functions and interpret the electron field as the wave function for one elec- 
tron in the lattice. In other words we would say that we are treating the problem in 
the mean field approximation where the operators are replaced by their mean values 
over configurations. It would be emphasized that this is not an exact calculation but 
just an artifact to obtain the best basis in which we would expand the operators of 
the Hamiltonian (2.1.1) in order to get the strong coupling regime. It can be easily 
shown that the following Lagrangean can generate the Hamiltonian (2.1.1) 



2\dt 2 ' V dt dt 
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_ ^drd± _ £^^*^y (2 1 2) 

2m dx dx a 

where ijj is normalized: 

j dx\^p{x,t)\^=l . (2.1.3) 
The equations of motion for the Lagrangean (2.1.2) are 

Equation (2.1.4) is the Schrodinger equation for an electron in a potential given by 

V(x, t) = — r](x, t) 
a 

while equation (2.1.5) is an equation for an oscillator with frequency oUo forced by 
the presence of an external field, | The picture is that of an electron which 

distorts the lattice which, in its turn, produces a potential for the electron; a self 
consistent interaction. 

We are interested only in stationary solutions for the electrons, that is, solutions 
of the form 

V'(a:,i) = (/),(x)e-^°*/^ (2.1.6) 
as well as in static solutions for the lattice (adiabatic solution) 

r){x,t)=r)o{x) (2.1.7) 
Prom (2.1.4), (2.1.5), (2.1.6) and (2.1.7) we get 
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We can think of (2.1.8) and (2.1.9) as follows: we put the electron in the lattice 
and the latter adjusts itself to the presence of the former. As consequence the 
electron is trapped by the potential well formed around itself. 

Substituting (2.1.9) in (2.1.8) we get 

- = E.^„ (2.1.10) 



2m dx'^ Mauj^" 

This is a non-linear Schrodinger equation which can be solved exactly. There is 
a localized static solution which reads 

^„M.i^,ec/.(^fcf^) (2.1.11) 



where 



and 

K--^ (2.1.13) 

8m 

is the binding energy of the electron in the potential well. Xo is an arbitrary constant 
which gives the center of the packet described by (2.1.11). 
The lattice displacements are given by (2.1.9) 

Vo{x) = -2a sech' ^9{x^^^ ^2.1.14) 

which is symmetric around the electron position. 

The potential well where the electron is trapped is given by: 

V(x) ^-2\Eo\ sech^ f 9{x-Xo) \ ^^.i.lS) 
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We can identify the parameter g as the strength of the interaction since the 
potential (2.1.15) becomes very weak for distances greater than g~^. Therefore, 
defines the polaron length. 

2.2 The Acoustical Case 

For electrons interacting with acoustical phonons the Hamiltonian can be written 
as i: 

- /"^a; I _j_ '^'"^ /^'^''^y _^ ^"^^ ^"^ _l_ £)^''^^t^ I (2 2 1) 

J I 2z/ 2 \dx ) 2m dx dx dx j 

where Vg is the sound velocity in the lattice and all other definitions are mantained. 

The second term in (2.2.1) comes from the Debye dispersion relation for acous- 
tical phonons 

uj = Vs\k\ (2.2.2) 

where k is the phonon wave vector. 

For (2.2.1) the "semiclassical" Lagrangean is 

h'^ dip* dip dr] 



D-^ri^} (2.2.3) 



2m dx dx dx 
which produces the following set of equations of motion 

■j^9ijj ^ %^ d'^ip ^ g (2 2 4) 

dt 2m dx"^ dx 

d'^T] ^d'^V D d \ ip ^ f2 2 51 

dt'^ dx'^ V dx 
These equations are interpreted in the same way as for the optical C3iSG clS cl self 
consistent interaction, where the potential that the electron feels is 
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Unlike (2.1.4) and (2.1.5), this new system of coupled equations admits a trav- 
elling solution. We can obtain a general solution for (2.2.4) and (2.2.5) defining a 
new variable, x = ^ ~ ^ "^t, where Xg and v are arbitrary constants. We look for 
solutions of the form 

ilj{x,t) = (j)o{x)exp \^^{mvx - E'j)^ (2.2.6) 

and 

v{x,t)=7io{x) (2.2.7) 
From (2.2.4), (2.2.5), (2.2.6) and (2.2.7): 

Eq. (2.2.8) is the Schrodingcr equation in the variable and (2.2.9) has the form 
of a wave equation of the same variable. Substituting (2.2.9) in (2.2.8) we get 

"'"'^^ ^'-^E'g-^]^o (2.2.10) 



2m dx^ V {vl -v'^y V " 2 
Again wc have obtained a nonlinear Schrodingcr equation which is characteristic 
of a self consistent interaction. The solution obtained if v is smaller than the velocity 
of the sound is 

UX) = ^sech (2.2.11) 

where 
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with the electron binding energy given by 
The solution is unstable for v > Vg. 

The interpretation is almost obvious: the electron and the lattice displacement 
move together with velocity v. Observe that the wave function (2.2.1) and the 
binding energy (2.2.13) have exactly the same shape as in the optical case for i> = 
(the adiabatic case) . We expect that the potential which the electron feels must be 
the same. From (2.2.9) 

and the potential is given by 

Observe that for i> = (2.2.15) has the same shape of (2.1.15). This exhibits some- 
thing fundamental in the physics of the two problems. Observe that the parameters 
g and g' has the same form for ^; = 

a \m) [WJ 

where Ec is the characteristic energy of the phonon system; hujo in the optical case 
and hcuo in the acoustical case, where 

i^D = Vsl a 

is the Debye frequency. So, g and g' play the same role in both problems and, from 
now on, we will call them g. Actually, we will work only with the static case, w = 0, 
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observing that this means that we are not taking in account the kinetic energy of 
the lattice in (2.2.3) (Born-Oppenheimer approximation). This explains why we are 
using the term "adiabatic" . 

Here, it should be noticed that one must reconcile g » 1 (strong coupling) with 
the continuum model we have been using. It has been shown in that there is a 
vast range of polaron sizes where both conditions are met. 

2.3 The Adiabatic Expansion 

As we have seen, for the acoustical as well as for the optical case, the potential 
well which traps the electron is the same. As a consequence the wave functions and 
binding energies are also the same. As we have treated the problem in a variational 
way we expect that the wave function, (2.1.11) or (2.2.11), is the ground state wave 
function for the electron in the adiabatic limit and zero temperature. Nevertheless, 
we expect that for finite temperatures the electron can be found in some excited 
state of this well due to its interaction with thermal phonons. As the interaction 
is self consistent the potential well must change its shape changing the potential 
energy of the electron. If the temperature is not too high so it does not remove the 
electron from the well, we might imagine that due to virtual transitions, the electron 
absorbs energy from the lattice and immediately emits this energy, remaining in the 
ground state. 

In order to find these excited states we have to solve the Schrodinger equation 
for the electron in the potential well, (2.1.15) or (2.2.15) (put Xo = 0) 




(2.3.1) 



This equation can be solved exactly and gives one bound state as expected 
(see (2.1.11) and (2.2.11)) 
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with energy 



E. 



8m 



and a set of doubly degenerate free states 



with energy 



where /c„ is solution of 



*„ = ^-^» = 0,±l,±2,... (2.3.4) 

and S{k) is the phase-shift due to the scattering of the free states given by 

S{k) = arctan [^^^) (2-3.5) 

Here we have imposed periodic boundary condition in x = ±L/2 with L ^ oo. 

The "adiabatic expansion" is made by expanding the electron field operators in 
(2.1.1) or (2.2.1) in this new basis, yielding 

CO 

i>{x, t) = (l>o{^)ao{t) + ^ 0^(x)a„(t) (2.3.6) 

n=l 

oo 

V^t(x, t) = MMit) + E ^n%^Kit) (2.3.7) 

n=l 

where a„ and aj^ are, respectively, the destruction and creation operators for each 
state of (2.3.1). They obey the following anticommutation rules 
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{an{t),am{t)} = {dl{t),al^{t)} = 

The field displacement as well as its conjugate momentum density can also be 
expanded in this basis as 



fi{x,t)=r]oix)+ Yl ^kit)^ (2.3.8) 

fc=_oo ViV 



k=—oo 
-oo gifex 



7r{x,t)= Pkit)^ (2-3.9) 



.=-oo v^^ 

where N — L/a is the number of ion sites and k is given by the periodic boundary 
conditions 

^ n = 0,±l,±2,... 



As f] and tt are real we must have 



Pk = P-k 



with the usual commutation rule 

[qk{t):Pk'{t)] = ih5k,k' 

If we substitute (2.3.6), (2.3.7), (2.3.8) and (2.3.9) in (2.1.1) and (2.2.1), and use 
the orthogonality of the adiabatic states 

I dx(f)o{x)(f)^{x) = (2.3.10) 

J dx4>l{x)4>l{x) = Snn' (2.3.11) 

and the equations (2.1.9), (2.2.9) and (2.3.1) we obtain 
i) for the optical case: 

oo 

eo + Eoo'a + ^ E^alan+ 

n=0 

12 



2^ + -^QkQ-k + -Qkif {k){a^a - 1)+ 



where 



k=—oo 



Y,{fn{k)a)an + /:(-A;)4a + ^ /nm(A;)at a^)]} (2.3.12) 

n=0 m=0 

(2.3.13) 

r(A;) = ^ / dx^l{x)e^'^^ (2.3.14) 
/„(^) = ^ / c?x0„(x)0^(x)e^^- (2.3.15) 
Um{k) = ^J dx<P^*{x)<P^{x)e^'^ (2.3.16) 



ii) for the acoustical case: 
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Ha = eo + EoO^a + ^ aj,a„+ 

n=o 

l^Jf + — ^ — IkQ^k + ikDqk[fo{k){a^a - 1) 

fc=— oo 

oo oo 

+ + fn{-k)<4fl + XI /nm(fc)a|iam)]} (2.3.17) 

n=0 m=0 

Notice that (2.3.12) and (2.3.17) are exact results, no approximations have been 
made so far. 

Since we are interested in the strong couphng hmit, the phonons can not excite 
the electron out of the well. Quantitatively this means that 



> 1 



where E^ is the phonon characteristic energy and | Eq \ the modulus of the bind- 
ing energy of the electron. The inequality above allows us to define an expansion 
parameter for the strong coupling limit, a, as 

Ec 

If now we scale the Hamiltonians (2.3.12) and (2.3.17) by the characteristic 
phonon energy we can easily see that the pure electronic part is of the order of 
while the interaction term is of the order of and consequently much smaller 

than the former 0, [0]. Using this result we will eliminate the electronic part of 
the Hamiltonian by a perturbative treatment and obtain a renormalized phonon 
Hamiltonian. 

Consider the pure electronic part as the non-perturbed Hamiltonian. Its eigen- 
states can be written in the Fock space as: 

I no,ni, ...noo > (2.3.18) 
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where Ho is the occupation number for the ground state. We simphfy the calculation 
assuming that there is only one electron in the problem (see appendix A for the case 
of many electrons), that is: 



j=0 



(2.3.19) 



The ground state is given by: 



V'o >=l = 1; rij = > j = 1, 2, 3, ... 



(2.3.20) 



with energy Eo. 

This state is exactly the adiabatic state given by (2.1.11) or (2.2.11). 

We can now apply the Rayleigh-Schrodinger perturbation theory to this problem. 
The first order term is null because wc have just one electron. The second order 
correction is easily calculated and it gives, 

i) for the optical case: 



D 



kk' 



where: 



n=0 



ii) for the acoustical case: 



Ef = -D^ kk'Voik, k')qm-k' 



kk' 



Therefore, we write the Hamiltonians for the renormalized phonons as 



(2.3.21) 



(2.3.22) 



(2.3.23) 



+ 00 



rr , V- ( PkP-k M ^ r o,A 



M-k'} (2.3.24) 



k=—oo 



k'=—oo 



where 
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e = eo + E, = (2.3.25) 

and 

[^1°,,)' = u'M, - ^Kik, k') (2.3.26) 

{nt)' = vlkH,,, - ^kk'Voik, k') (2.3.27) 

The Hamiltonian (2.3.24) was already obtained some years ago for the acous- 
tical case and it can be rewritten in a different form in terms of some new operators 
6fi{x,t) and 7r{x,t) defined by (c.f. (2.3.8)) 

q,(t) = ^ f -e-^'^^Sr^ix,t) (2.3.28) 



Pk{t) = I dxe-''''^7f{x,t) (2.3.29) 



/N 

One can easily check that 

[Sfi{x, t), 7r(x', t)] = ih6{x — x ) 
Now, using (2.3.28) and (2.3.29) in (2.3.24) we get 

{"2 2 2 

|- + - ^5fj j dx'F{x, x')5Ti{x') (2.3.30) 

and 



where 

F{x, x') = 4g M^)M^') E fjiy^)^!^ (2.3.32) 
These results agree perfectly with those obtained by other methods 0, 
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3 The Effective Hamiltonian 



As we are interested in the excitation spectrum of the Hamiltonians (2.3.30) and 
(2.3.31), we have to diagonahze them. If we choose 



5ri{x,t) = J2 qn{t)un{x), (3.1) 

n=0 

where Un are the normahzed eigenfunctions of (2.3.30) or (2.3.31) 

r dx 
J a 

they must satisfy the following integro-differential equations 



U*^{x)Um{x) = Snrn, (3.2) 

a 



l) Vtlun{x) = UjlUn{x) - j dx' F {x, x')Un{x') (3.3) 

for the optical case, while 

^^) ^>r.{x) = -vl'^^ + vljdx'^F{x,x'f-^ (3.4) 

for the acoustical case. 

Observe that to solve (3.3) and (3.4) is equivalent in the momentum space to 
diagonalize (2.3.26) and (2.3.27). 

These equations often appear in works about polaron dynamics [0], 0, and 
we do not intend to solve them in this article. For (3.3) there is a closed solution 
[P] which we will use later, while for (3.4) we will make some approximations which 
are suitable for our purposes. 

An interesting and important solution of those equations can be found directly. 
These are the zero mode solutions, that is, solutions with fi^ = (see appendix B). 
For the optical case we have 



w^^(x) = ^^^tanh ( seek' ( ^ ) (3.5) 
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and for the acoustical case 



ufi.) = ^.ecA= (f ) (3.6) 
Examining (3.5) and (3.6) and comparing with (2.1.14) and (2.2.14) we see that, 



in both cases, 



Uo{x) = C ^Vo{x) (3.7) 

where C is a constant which appears due to the normahzation of Ug- The above 
relation clearly expresses the translational invariance of Hamiltonians (2.1.1) and 
(2.2.1). Notice that although we have put = in (2.3.1), in order to obtain the 
adiabatic basis, the center of the soliton solutions (2.1.11), (2.1.14), (2.2.11) and 
(2.2.14), is arbitrary and therefore we must have the m„ also expanded about this 
point, that is: 

Un = Unix - Xo) (3.8) 

Suppose that we move the center of the functions (2.1.14) or (2.2.14) by an 
infinitesimal quantity 6xo-Then, 

U ix jOX 

Tjoixo + Sxo) ~ rioixo) + -^-^ — - (3.9) 

where we have used (3.7). By (3.9) we conclude that the zero mode frequency 
corresponds to the translation of the soliton, in other words, to the motion of the 
polaron. 

Once we have the eigenfunctions of (3.3) and (3.4), we would expect to write the 
Hamiltonian in the form: 



n=l 



where we have used that fio = 0. 
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At first sight (3.10) shows a free particle with momentum Po and a set of de- 
coupled harmonic oscillators. Nevertheless, it is not possible to take it seriously 
because, initially, we have implicitly assumed that the lattice displacement can not 
be indefinitely large. From (2.3.21) and (2.3.23) we see that the energy correction 
depends on the lattice displaceme nt which must be finite in order to the pertubation 
theory to be valid. 

Let us observe that, due to (3.9), the polaron displacement is proportional to 
the displacement of its center, that is: 

qo = 6xo/C (3.11) 

From (3.1), we have 

oo 

5ri{x, t) = QoUo + QnUn 
n=l 

or using (3.7) and (3.11) 

6r]{x,t) = 6xo-^ + 'YqnUn (3-12) 

Therefore we can assumed that true dynamic variable, that is, Xo = Xo{t). 

So, based on (3.12), we will rewrite the expansion (3.1) as: 

oo 

(5r/(x, t) = r]o{x - Xo{t)) + qn{t)un{x - Xo{t)) (3.13) 

n=l 

This procedure is known as "collective coordinate formalism" |1^ . Observe that 
(3.13) changes the kinetic part of (2.3.30) or (2.3.31) because Xo is also a function 
of time. It is shown 0| that the new Hamiltonian in the presence of the polaron 
position operator, Xo, is given by 
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8Mo 



nniQnPn 



n,m=l 



n=l 



(3.14) 



here {, } denotes anticommutation and P is the momentum operator associate to Xg 



and 



[xo{t),P{t)] = ih 



[qn{t),Pm{t)] = ihd„ 



with all the other commutators being zero. 
In (3.14) we have 



drio 



dx 



/+00 
dx 
-oo 

as being the classical soliton mass which becomes 



Mr, 



8 KflUJr 



for the optical case, and 



M„ = — m 



3 Vhvsg ^ 



for the acoustical case. 
The new quantities 



M r dx drjo du.. 



Mo J a dx dx 



(3.15) 



(3.16) 
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and 



G'_ = /^«™(^)^^ (3-17) 

couples the polaron to the renormahzed phonons. 

In the strong couphng hmit (3.14) can be simphfied f^. Due to (3.5), (3.6), 
(2.1.14), (2.2.14) and (3.7), we can rewrite Sn as 

Ml fEo\ f dx dUn, . 



Mo y/ag \ D J J a dx 
As the integral only gives a numerical factor, this yields 

M 1 fEo\ 1 



Sn 



Mo ^/ag \ D J a 



Now, from (3.15), 



and 



M M /E^^ 



{ME,Y/2 



so 



Snln 



Er 



.1/2 



< 1 



Eo 

and therefore, this sum is very small in the strong coupling limit. Within this 
approximation we get 



* - + E (ij + -^€) + - ^ (3.18) 

The second term in (3.18) is the energy of non-interacting phonons and the third 
term can be interpreted as the kinetic energy of the polaron. Observe that 
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1 - 1 - 

•^o — "TIT [2^0 -^^1 — , r I P / , G nm^.mVn 
in Mo ' 



(3.19) 



n,m=l 



and so, P can not be the polaron momentum because since 

it is a constant of motion. Prom (3.19) we interpret M^Xo as the polaron momentum 
and GnmQmPn the momentum of the phonon field. Observe that Hamiltonian 
(3.18) is very close to the electromagnetic Hamiltonian, where the coupling between 
the particle and the field are obtained via the potential vector (see Section 4). 
If we define the destruction and creation operators 



V 2h 



MQr, 



2h 



which obviously obey 



(3.20) 
(3.21) 



[K, bl 



mi '-'nm 



one can rewrite H as 



where 



00 1 
H^e+Y: hn„{b% + 1/2) + +^{P - Pphf 



n=l 



ph 



E - . 



00 *^ 

+ y - 



(3.22) 



Here we have used the fact that, from (3.15), is anti-symmetric in the inter- 
change of m and n. 
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Gnm — —Gjnn (3.24) 

Observe that the momentum of the phonon field consists on two parts; a diagonal 
part (the first term on the right hand side of (3.23)) which commutes with the 
phonon-number operator: 

oo 
n=l 

and, therefore, conserves the number of phonons in the system. This term is respon- 
sible for scattering. The other term does not commute with the number operator 
and is related with absorption or emission of phonons by the polaron (Cerenkov pro- 
cess) . We will restrict our problem to typical polaron kinetic energies much smaller 
than the phonon energies, in other words, small velocities. In this limit the occur- 
rence of emission or absorption of phonons is not possible due to the simultaneous 
conservation of momentum and energy. In terms of our parameters this means that: 



(3.25) 



Only scattering, and therefore virtual transitions, will be relevant for our problem. 

With this approximation the polaron dynamics will be described by the following 
effective Hamiltonian: 



2 

oo 



^ = ^ I ^ - E ^9nJlk 1 + E ^^nhik (3.26) 



\ n,m=l / n=l 



where 



ynm — „. /rwT — ^nm K'-'-'^') 

As we will show in the next section the Hamiltonian (3.26) describes the dynamics 
of a brownian particle, that is, a heavy particle in a bath of light particles which 
collide with it. 



23 



4 Functional Integral Method 



The starting point for the calculations of the transport properties of the polaron 
is the well-known Feynman- Vernon formalism that the authors have recently 
applied to the Hamiltonian (3.26). 



We are interested only in the quantum statistical properties of the polaron and 
the phonons act only as a source of relaxation and diffusion processes. Consider the 
density operator for the system polaron plus phonons, p{t). This operator evolves 
in time according to 

p(t) = e-'^'/^p{0)e'^'/^ (4.1) 

where H is given by (3.26) and p(0) is the density operator at t = which we will 
assume to be decoupled as a product of the polaron density operator, ps{0), and 
the phonon density operator, p_r(0): 

p(0) = ps(0)p^(0) (4.2) 

where the symbol S_ refers to the polaron (system of interest) and R to the phonons 
(the reservoir of excitations). 

The condition (4.2) means that we put the electron in the lattice which is in 
thermal equilibrium at temperature T. So, we consider the phonons as described by 
their equilibrium distribution, 

Pr{0) = (4.3) 

where 

Z = tTR (e-^^^) (4.4) 

with 
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Here trR denotes the trace over the phonons variables and Kb is the Boltzmann 
constant. Hb is the free phonon Hamiltonian which is given by the last term on the 
right hand side of (3.26). 

As we said, we are interested only in the quantum dynamics of the system S, so, 
we define a reduced density operator 

Ps{t)=trn{p{t)) (4.6) 

which contains all the information about S when it is in contact with R. 

Projecting now (4.6) in the coordinate representation of the polaron system 

Xo I Q >= 1 \ Q > (4.7) 
and in the coherent state representation for bosons (the phonons) 

bn I a„ >= a„ I a„ > (4.8) 
we have |]l3l (see also appendix C) 

Ps{x, y^t) = J dx' j dy'J (x, t] x\ y' , Q)ps{,x' , y\ 0) (4.9) 

Here we have used (4.1), (4. 2), (4. 6) and the completeness relation for the repre- 
sentations above, namely, 

J dq\q><q\=l (4.10) 

f d^a , , , , 

/ |a><a|=l (4.11) 

where d'^a = d{Rea)d{Ima) as usual. 

In (4.9), J is the superpropagator of the polaron, which can be written as 
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J- Tdx e^(^°l^l-^°[^l)F[x,y] (4.12) 

Jx' Jy' 

where 

SM-jj^[^} (4.13) 
is the classical action for the free particle. F is the so-called influence functional 

^, , f (fa f (f(3 f (fp' _ic;|2 _i£!i_iZi! 

/"* D^a r D^l e^^[^''^]+^?[^'^] (4.14) 

where (3 denotes the vector {Pi, (32, Ps, ■■■,Pn) and Si is a complex action related to 
the reservoir plus interaction 



S,l.. S] ^ { *' |i (a . f - S- . f ) - - ih,) } (4.15) 

with 



HR=J2^^na::,an (4.16) 

n=l 

oo 

^/ = X) %gnmOi*man (4.17) 
n,m=l 

Here we have obtained a result which is very close to the electromagnetic coupling 
where the Hamiltonian depends on the vector potential, A, through 

^\ 2 

eA 

\P 

V ^ 

but the Lagrangean depends on 

V ■ A 
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In our case the Lagrangean formulation simplifies the problem transforming a 
nonlinear problem into a linear one. The action (4.15) is quadratic in a, so it can 
be solved exactly. Observe that the Euler-Lagrange equations for (4.15) are 



din + i^nOCn " X] 9mnOlm = (4.18) 
m=l 
oo 

d* - iO^nOin + ixY^ 9nmOCn = (4.19) 

m=l 

which must be solved subject to the boundary conditions 

«„(0) = /?„ (4.20) 

<it) = < (4.21) 

Due to (3.24) we have = 0, so, the modes are not coupled among themselves. 
This makes (4.18) and (4.19) easy to solve. That set of equations represents a set 
of harmonic oscillators forced by the presence of the polaron. The result can be 
written as: 



«n(r) = e-^^"" /3„ + E Wnm{T)Pm (4.22) 
\ m=l / 

:(t) = e^""^ f<e-'^"* + f: Wnm{r)e-'''-'aA (4.23) 



where Wnm ^-nd Wnm s-re functionals of x{t) which obey the following equations 

oo 

Wnm - + E W:,,Wn'm (4.24) 

n'=l 

oo 

Wnm = + E W^n'Wn'm (4.25) 



n'=l 



where 



W:M,r) = t r dt'gnmHty^''''-''-^'' (4.26) 
Jo 
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r) = i [ dt'grnr.xity^''--''-^'' (4.27) 

(observe that: Wnm{t) = W^mn(O)). 

Now we expand the action (4.15) around the classical solution (4.22) and (4.23) 
and obtain, after some integrations in (4.14) 

CO 

F[x, = n (1 - r„„[x, y] rr„)-^ (4.28) 

n=l 

where 

oo 

^nm = W*Jy] + Wrnnix] + E ^/mbl^^^nW (4.29) 

with 

n,^^ = (e'^^^" - l)~^ (4.30) 
Notice that (4.28) and (4.29) are exact, no approximations have been made so 

far. 

We see from (4.24) and (4.25) that can be expressed as a power series of the 
Fourier transform of the polaron velocity, i, so, due to the small polaron velocity 
condition (3.25), we expect that only few terms in (4.24) will be sufficient for a good 
description of the polaron dynamics. 

Another way to sec this, is to notice that (4.24) and (4.25) are the scattering 
amplitudes from the mode k to the mode j. The terms that appear in the sum rep- 
resent the virtual transitions between these two modes. With these two arguments 
in mind we will make use of the Born- approximation. In matrix notation, 

W ^{l-Wy^W + W°W° (4.31) 

Therefore, in the approximation of small polaron velocity the terms in (4.29) are 
small and we can rewrite as a good approximation 
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F[x, y] ~ exp | ^nn[x, y] rr„| (4.32) 

Observe that if the interaction is turned off (F 0) or the temperature is zero 
(T = 0) the functional (4.29) is one, and, as we would expect the polaron moves as 
a free particle. 

Substituting the Born approximation (4.31) in (4.32) and the latter in (4.12) we 
find 



J- £dx l''Dyexpl^^S[x,y] + ^^[x,y]^ (4.33) 



where 



S^l dt'[-^{x\t')-f{t')) + {x{t')-y{t')) I dt" Tj{t'-t"){x{r)+y{t"))'j 

(4.34) 



and 



^= fdt' fdr{rn{t'-r){x{t')-y{t')){x{r)-y{r))} (4.35) 
Jo Jo 



with 



TR{t) = hOit) £ gl^uncosiyt^ - njt (4.36) 



n,m=l 

oo 



ri{t) = he{t) glmnnsin{n^-nm)t (4.37) 

n,m=l 

Now, if we define the new variables R and r as 

R = ^y (4.38) 

r^x-y (4.39) 
the equations of motion for the action in (4.34) read 
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R{t) + 2 /* rft'7(r - t')R{t') = (4.40) 
Jo 



f (r) - 2 / dt'-f{t' - r)f(t') = (4.41) 



where 

or, using (4.37) 

7(t) = E ^7L^n(^^n - ftm)cOs{Qn " ^^«^)^ (4.43) 

^'^o n,m=l 

is the damping function. 

In terms of these newly defined variables, we can easily see that (4.40) and (4.41) 
have the same form of the equations previously obtained in the case of quantum 
brownian motion |T^, except for the fact that they now present memory effects. It 



should be emphasized that although (4.40) and (4.41) have only indirect physical 
meaning, through the study of the motion of the center of a wavepacket and the 
spreading of its width, 7(t) really plays the role of the damping parameter in the 
equation of motion of the former (see ref. [|1^ for details). 

Furthermore, we shall prove that (4.43) can be written in the form 



7(t)=7(T)5(t) (4.44) 

where 7(T) is a damping parameter which is temperature dependent and 6{t) is the 
Dirac delta function. The form (4.44) is known as the Markovian approximation, 
because in this case the memory is purely local and does not depend on the previous 
motion of the particle. 

If we use (4.40) and (4.41) with (4.44) and expand the phase of (4.33) around this 
classical solution we get the well-known result for the quantum brownian motion [|I^ 
where the damping parameter 7 (temperature independent) is replaced by 7(T) and 
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the diffusive part is replaced by (4.35). As a consequence, the diffusion parameter 
in momentum space will be given by 



D{t) = h—f = -n^e{t) J2 alnM^n - nmfcos{nn - nm)t (4.45) 

n,m=l 

We will also prove that D{t) has the Markovian form: 

D{t) = D{T)6{t) (4.46) 

where D(T) and 7(T) obey the classical fluctuation-dissipation theorem at low tem- 
peratures |1T5[| . 

In what follows we shall define a function S{uj,uj') which will, in analogy to 
the spectral function J{uj) of the standard model []T^, allows to replace all the 
summations over k by integrals over frequencies: 

oo 

S{u, J) = 9lJ{^ - ^n)5{uj' - n^) (4.47) 

n,m=l 



Notice, however, that unlike J (to) in this new function S{uj,uj') is related to 



the scattering of the environmental excitations between states of frequencies uj and 
u' (as seeing from the laboratory frame). Morover, due to (3.24) it is easy to see 
that 

S{uj,uj') = S{uj',uj) (4.48) 

From now onwards we shall call S{uj,uj') the "scattering function". 
Notice that we can rewrite (4.43) and (4.45) as 
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TiO[t) [°° , , 

^(^) = '^TT / / dLj'S{uJ,u'){uJ - u'){n{uj) - n{uj')) 
ZMo Jo Jo 

cos{u - uj')t (4.49) 

and 

2 Jo Jo 

{n{uj) + n{uj'))cos{uj - uj')t (4.50) 

Concluding, we have established that the Hamiltonian (3.26) leads to a brownian 
dynamics, that is, the polaron moves as a particle in a viscous environment where 
its relaxation and diffusion are due to the scattering of phonons. 

5 MOBILITY AND DIFFUSION 

Equations (4.43) and (4.45) show that the polaron transport properties depend 
essentially on the coupling parameter Qnm- From (3.17) we see that this parameter 
can be obtained if we know the eigenfunctions of (3.3) and (3.4). 

First of all we can show that (3.3) and (3.4) have solutions with definite parities. 
This is easily seen by changing x by —x in (3.3) and (3.4) and x' by —x' in the 
integral term. From (2.3.32) we observe that F{ ') and therefore 

Un{x) and «„(— x) obey the same eigenvalue equation. In other words, the Hamilto- 
nians commute with the parity operator and therefore it is possible to classify their 
eigenfunctions as odd or even. Now we must study the optical and the acoustical 
cases separately. 

5.1 Optical Case 

Turkevich and Holstein obtained the exact solutions for (3.3). For the odd modes 
the eigenfunctions are: 
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n = 0,2,4,6... 



where 



Y{x) = tank (5.1.2) 

and Pn are the Legendre polynomials. 
The eigenvalues of the problem are 

In particular the zero mode, n — and fio = 0, is given by (3.5). 
The even modes can be written as 



1/2 

, . fag f 2a + 5 V (1-^^(3;)) d 
"°<") = VTi (a + 2)(a + 3) j (5.1.4) 

where the allowed values of the a are solutions of 



V^(a + 3) - = -tan ^—j 



and ip is the digamma function. Its eigenvalues are given by 

^"^"4^- cv^ + 5a + 4 j 

We will use the convention given in table 1. 

As the labels for the even solutions are not integers we define, e„ — n — a, as 
the difference between our classification and the label. Table 1 shows to us that the 
eigenfrequencies go quickly to cUg while e„ goes to zero. 

Prom (3.17) we note that Gnm only couples functions with opposite parity. Sub- 
stituting (5.1.1) and (5.1.4) in (3.17) we get: 
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Tl 








n 

yj 


n 


n 


n 


1 


0.523 


0.648 


0.477 


2 


2 


0.882 





3 


2.601 


0.912 


0.394 


4 


4 


0.949 





5 


4.648 


0.958 


0.352 


6 


6 


0.971 





7 


6.674 


0.975 


0.326 


8 


8 


0.981 





9 


8.692 


0.983 


0.308 


10 


10 


0.987 






Table 1: Conventions for classification of the eigenfunctions 



^ _ 2^ gm(7r(n - m + e^)) 
TT [(n - m + e„i)^ - IJ 



where 



_ [(n + 2)(n + 3)(2n + 5)(m - + 2)(m - + 3)(2m - 2e„ + 5)]i/2 



[{n + m + emf + 10(n + m - e^) + 24] 

(5.1.6) 

for: n = 0,2,4,6... m = 1,3,5,7... 

Observe that (5.1.5) is strongly peaked around n = m ± 1. Therefore the most 
important contributions to summations involving Gnm will come from these forms 
(observe that goes to zero as m goes to infinite). 

Gnm - gKnm[S{n - 171 - 1) - S{n - III + 1)] (5.1.7) 

for n even and m odd. 

From (4.47) and (3.27) we get: 

S{u;,u;') = -^^{C^„_J5(a; - n„)5(a;' - Qn-i) + 
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+ (5.1.8) 

where n is even and 



Substituting (5.1.8) in (4.48) we find: 

^ "slr^^^) E{Cnn-l(^n - ll„-l)(n(Q„) - n(J^„_i))cOs(Q„ - nn-l)t 

+ Ci+i(Q„+i - Q„)(n(Q„+i - n(Qj)cos((Q„+i - Qjt)} (5.1.10) 
We will define wavevectors for each n in (5.1.2) in the form 



= f (5.1.11) 

where L is the length of quantization (L — > oo). (5.1.10) then becomes 



^(t) = Ji^e(t)— / dKC^iK, K)\ — \ Utt cos -^t 5.1.12 

^ 8Mo ^ ' L Jo ^ ' ' \dK \dQ \LdK ^ ' 



where we used the limit L — > oo. 
It is easy to see that by (5.1.3) 



dK L2 
and by (5.1.9) and (5.1.6) that 



dVt A-k"^ ijJo 

as L ^ oo (5.1.13) 



r2 1^2 

C\K, K) = — — asL^oo. 



7r2 



Now, using the fact that the frequencies approach very fast the value ujq when n 
increases we make the following approximation: 
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dn dn 

we shall rewrite (5.1.12) as 

This integral can be easily done if we change variables, x — The integral 

then becomes 

and finally 

'^^ 2Mo3KbT V y ^ ^ 

which has the form (4.44). 

For low temperatures, KbT <^ hwo, we have 

so 7(7") goes to zero as T goes to zero, as expected. So, for very low temperatures 
the mobility is extremely high and the polaron moves as a free particle. This is an 
expected result since at T = there are no phonons to be scattered. 
For high temperatures, KbT ^ hwo, 



6 Mo V hojo 



and the mobility goes to zero as T — > 00. 

The diffusion parameter (4.49) can be calculated in the same way. It gives 



^(r) = ^^^a;„(e^°/^-^-l) (5.1.17) 
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For low temperatures this parameter is too small, going to zero as T — 0. 
The fluctuations are once again small due to the absence of phonons. For high 
temperature we see that the fluctuations also increase linearly with T. 

Observe that 



liT) 



2MoKbT (l - e-^'^°/^sT^ (5.1.18) 



which gives the classical result of the fluctuation-dissipation theorem [T^ for the 
brownian motion at low temperatures. 

5.2 The Acoustical Case 

As we do not have exact solutions for (3.4) it will be necessary to make some 
approximations in the present analysis. Observe that (3.4) is a Schrodinger-like 
equation for a particle in a non-local potential 

d'^F{x,x') 



V{x, x) 



dx'dx 



From (2.3.32) we observe that V{x, x') goes to zero as x goes to infinity. Actually, 
the potential is almost zero except in the range 



1 1 

- < X < - 
9 9 



Out of this range the wave function can be well described by 



dx"^ 



+ rMfc(x) = (5.2.1) 



where we have used that 



u: = Vs\k\ (5.2.2) 
The solutions of (5.2.1) must be classified as even or odd. We choose 
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UE{k,x) = \ —cos{k I X I +SE{k)) 



uo{k,x) 



where 



sgn[x) 



sgn{x) sin{k \ x \ +5o{k)) 



1 if X > 
-1 if X <Q 



(5.2.3) 
(5.2.4) 



and Ssik) and So{k) are the phase shifts for the even and odd modes, respectively, 
which must appear due to the presence of the potential 
Another possible solution of (5.2.1) is: 



) = [t{k)e'''''e{x - 1/g) + (e*^^ + r(A;)e-*^^) e{~x - 1/g)} (5.2.5) 



Uk(X 



This expression can be interpreted as a wave incident from the left on a potential 
whose t{k) and r{k) are the transmition and reflection amplitudes. 
We can construct (5.2.5) from (5.2.3) and (5.2.4) as |jl6| 



Uk{x) = e'^''uE{k,x) +ie'^°uo{k,x) 



if 



m = I ( 



^2iSE{k) _^ ^2i5o{k) 



1 



Consequently the transmission and reflection coefficients are given by 



T{k) =1 t{k) |2= cos\6oik) - 6E{k)) 
R{k) =1 r{k) p= sin^{5o{k) - 6E{k)) 



(5.2.6) 
(5.2.7) 
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and T{k) + R{k) = 1 as expected. 

Once we have the phase shifts of the problem we can find the transmission and 
reflection probabilities using (5.2.6) and (5.2.7), or alternatively, if we have the 
reflection and transmission amplitudes we can obtain the phase shifts 



6E{k) 



1 



-arctan 



6o{k) = -arctan 



Im{t{k) + r{k)) 
Re{t{k) + r{k)) 

Im{t{k) - r{k)) 



(5.2.8) 
(5.2.9) 



Re(t{k) - rik)) 

Actually, Schiittler and Holstein |0] obtained these coefficients in the limit of 
long and short wavelength after a rather intricate algebra: Q 
i) ioi k^ g 



r{k) 



ii) for k <^ g 



r{k) 



3ik 


ge 


9 


9' 


3ik 


9k' 


9 


9' 



These results allow one to compute the respective phase shifts as: 
i) for k ^ g 



(5.2.10) 
(5.2.11) 



(5.2.12) 
(5.2.13) 



SEik)c^Soik) 



9_ 

5k 



ii) for k <^ g 



(5.2.14) 



^ In the results of ref. 5, R and T depend on the polaron velocity which is very small in our 
case and we have put it equal to zero. 
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So{k) ~ (5.2.15) 

So, the phase shifts are very small. We would say that there is a propagation of 
sound waves through the polaron. 

If the interaction between the electron and the lattice is strong, the range of the 
potential is small. Therefore, the contribution to the integral in (3.17) due to the 
true solution is almost the same as the one we would have got had we used the free 
solutions (5.2.3) and (5.2.4). 

Firstly, we impose periodic boundary conditions which gives the allowed values 
for k: 

2n7i 

kn = —j- n = ±1, ±2, ±3 . . . (5.2.16) 
In order to classify the solutions we will use the following convention: 

U2n-l{x) = UE{n,x) (5.2.17) 

U2n{x) = uo{n,x) (5.2.18) 

for n = ±l,±2,±3... 

Now we can evaluate (3.17). It yields 

2kn fSin[{kn ^ km)L/2\ /c /, \ r n w 
G2n-i,2m = ^{ , _ , cos(5e(A;„) - do{km)) 



+ 



1 - cos[{kn - km)L/2] 



sin{SE{kn) - So{km))} (5.2.19) 



k — k 

As in the optical case we have a matrix with zeros in the diagonal and with off- 
diagonal terms which decrease as a function of their distance to the main diagonal. 
When L ^ oo we will have (using (5.2.6)) 
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(5.2.20) 



where P denotes the principal value. 

Substituting (5.2.20) in (4.47), transforming the summations into integrals and 
using (5.2.2) we get 

SM = ~^^.^^)S{u;-u;'y^^^!^^^ (5.2.21) 

In (4.48) we will change the variables of integration and rewrite (4.47) as 

n0(t) r^D fi^D / Q Q 



^ 2MoJo J-u^u V 2' 2 7 

^ §) ~ ^ ~ §)) (5.2.22) 



where 



Q. — uj — uj' 
2 

Observe that we have replaced the limit on the integration by the cut-off fre- 
quency, UJd- 

Actually we are interested in a time scale, r, which is much longer than the 
typical phonon period or 

With this approximation the cosine term in (5.2.22) oscillates rapidly, giving 
no contribution to the integration, except when VL is close to zero. So we can 
approximate (5.2.22) as 
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where 



m = -hm/s{9+'-,e-'-^ (5.2.24) 

Now, using (5.2.21) and noticing that the delta term does not contribute to (4.48) 
we get 

= T^^'^W (5-2-25) 

where we have used (5.2.7). 

So, we conclude that we have here a Markovian process with the damping pa- 
rameter given by 

7(T) = , ^ / duu'R{u) (5.2.26) 

Defining a new variable 

uo = gVsK,/2 
and a typical phonon temperature, Tc, by 



2Kb 

we can evaluate (5.2.26) which reads 



327rMo V T 

where 



nS) = Sj^ dK R{k) _ (5.2.28 



is exactly the result obtained by Schiittler and Holstein for the polaron mobility 
using the kinetic theory. 
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For small temperatures, T <^ Tc, we shall use the long wavelength reflectivity 
(see (5.2.12)) since it gives the largest contribution to the occupation number 



R{k) ~ (5.2.29) 
Then, using (5.2.28) and (5.2.29) we can approximate (5.2.25) by 

This result shows that the acoustical polaron, as the optical one, behaves as a 
free particle as T ^ 0. 

For high temperatures, T » we use the expression for short wavelengths (see 
(5.2.10)) 

R{k) ~ ATT^K'^e-^'"' (5.2.31) 

and one has 

which means that the mobility decreases for high temperatures. 

We can calculate the diffusion coefficient (4.49) in the same way and we get 

where 

and we have used the fact that the diffusion is a Markovian process. 
For small temperatures, T -C the diffusion coefficient is given by: 

D{T) ~ ^I^Kb'^ (5.2.35) 
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and the fluctuations decrease very fast as tlie temperature is lowered, exactly as in 
the optical case. 

So, the relation between relaxation and diffusion is the classical one for the 



brownian motion 15 



7(T) 

For high temperature, T ^ Tc, 



MoKbT 



(5.2.36) 



(5.2.37) 



And exactly as in the optical case the fluctuations increase linearly with tem- 



perature, this is the classical result |jT5| which is expected to be valid in the high 
temperature limit. 



6 CONCLUSIONS 

In the foregoing sections we have shown that the semiclassical (mean field) method 
enables us to visualize the polaron physics and allows us to treat the strong coupling 
limit of an electron interacting with a lattice. 

The advantage of dealing with this method is the fact that, in terms of the 
coordinate and the modified phonons, we reduce the problem to a new model for 
treating quantum dissipation. In a sense, the non-linear character of the electron- 
phonon interaction is somehow "hidden" in the soliton-like solution whose center is 
regarded as the polaron coordinate. 

Eliminating the electron operators by perturbative techniques (that is, tracing 
over the electron coordinates ) and using the well known collective coordinate formal- 
ism we get an effective Hamiltonian for the polaron in the presence of renormalized 
phonons. That Hamiltonian, in the approximation of small polaron velocity, is re- 
duced to a very simple form which takes into account only processes which involve 
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polaron-phonon collisions. 

We developed a new functional method to treat the Hamiltonian in the limit of 
small polaron's velocity. Our method showed that the polaron moves as a brownian 
particle which collides with the light particles of the environment. This method pro- 
vided us with a tool for a systematic calculation of the damping parameter (and, as 
a consequence, the mobility) and the diffusion coefficient as function of the tempera- 
ture. We have also showed that in the time scale of interest the motion is essentially 
Markovian, that is, it does not have memory. 

An important comment about our work is that it is fully quantized and the 
"semiclassical" argument is only used as an artifact. Furthermore, it confirms some 
important results for the acoustical polaron obtained by Schiittler and Holstein 
using kinetic transport theory. 
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de Amparo a Pesquisa do Estado de Sao Paulo) for a scholarship, while (A.O.C.) 
kindly acknowledges the partial support of CNPq (Conselho Nacional de Desenvolvi- 
mento Cientifico e Tecnologico) and FAEP (Fundo de Apoio ao Ensino e Pesquisa 
da UNICAMP). 
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APPENDIX A 



In this appendix we wish to calculate the first order correction in energy due to a 
many electron wave function in (2.3.12) and (2.3.17). 

In case of many electrons the ground state is the Fermi sphere with radius kp, 
which is given by 



2a \N 

where is the number of electrons and N the number of sites. 
The non-perturbed Hamiltonian is 

Ho^ Eoa^a + J^Enaian (2) 

k 

with the ground state wave function 

I ip° >=\ Uk = I, k < kp; Uk = 0, k > kp > ■ (3) 
The ground state energy is 

where we have accounted for the spin degeneracy. 

The first order correction is given from the interaction term in (2.3.12) or (2.3.17) 

E^'^ =< ilJo \ Hi \ ilJo > 

For the optical case it reads 

E^'^^-J2f,,{k)e{kp-\q\)qK 

" qk 

or 
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^(,)^ A(fc)^^ (5) 
T « 

where 

For the acoustical case, 



E^'^^tJ2kA{k)qk (7) 
So, at first order the ions are displaced from their equilibrium positions by 



A{-k) 

in the optical case and 



Macul) 



ikA{-k) 
MujI 

in the acoustical case. This causes a change in the energy given by 

A(fc)A(-fc) 

for the optical case, and 
for the acoustical case. 

APPENDIX B 

Let us first show how to obtain the zero mode solution for the optical case. Making 
Qo = in (3.3) we get 
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uo{x) = j dx'F{x,x')uo{x') 



Define a function g{x, x') 



g{x,x) = 2^ 

n= 



Then, from (2.3.1) 



= -^sech (^y) ^^c^ (^^) + '^(^' 
where we have used the completeness of the adiabatic states 



n=0 

and the exphcit form for 4>o{x)- 
Let us rewrite (1) as 



where 



uo{x) = g'^sech (^y^ /(x) 



Then, from (3) 



/(x) = dx'g{x,x')sech ( ^ ) -uola^') 



5(X 



= sech —Mo (2^) + ^ dx'sech 

The last term in (7) must vanish because 1*0(2^) is odd. Using (5) we 



48 



The solution is easily obtained and reads 



/(x) = sech tanh 
Substituting (9) in (5) and normalizing it, we get 



, , / 15ag f gx\ , 2 f 9^ 
uo{x) = W -^tanh \ — 1 sech [ — 

For the acoustical case we must use (3.4) with VLq = 0, which reads 

(fuo d f , , duc 



dx^ dx 

Now, using (2) we should define 



/ dx'F{x,x')'^{x') 
J dx' 



h{x) = J dx'g{x,x')sech ^^^^ "(f^^^'^ 
and rewrite (11) as 

d^uo 2 d r ,f9x.,. ^ 
This can be easily integrated yielding 



dx' sech ( j 



h{x) 



where 



Mo(— oo) = 



Observe that h{x) also obey eq. (8), so 



h{x) = tanh ^^ch 
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Substituting (14) in (13) we get, after normalization 



'^o{x) = y —sech^ (^yj (15) 

APPENDIX C 

We shall evaluate here the functional form for the superpropagator, J. Prom (3.26) 
we see that the Hamiltonian can be put in the form 

H^Hs + Hr + Hj (1) 

where 

Hs = Py2Mo (2) 

oo oo 

Hr = J2 fi^nbibn + ( E figr,mhlhr,f/2Mo (3) 

n=l n,m=l 

oo 

Hi = -P E hgnmblbjMo (4) 

n,m=l 

Prom (4.6) we get 

< X I Ps(*) \y>^ J J (n < ^« I pit) I y« > 

or using the completeness relations (4.10) and (4.11) we can, with the help of (4.2), 
write 

Psix,y,t) = j dx' j dy' ps{x',y',0)j{x,y,t;x',y',0) (5) 

where 

d'^a fd'^0 fd^p' 



-^ = 1^1^! ^P^^^*^ ^Wxa*x'P, t)K*{yayf*t) (6) 



with 
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K{xa*x'(5t) =< xa \ e''"*'^ \ x' (5 > (7) 

In order to transform (7) into a functional integral we must divide t in (M-1) subin- 
tervals of length e and use (M-1) completeness relations between the (M-1) expo- 
nentials in (7). Then 

< xa I e-^*/^ I x'(5 >== J dq.., . . . J dq, J ^ . . . / ^ 

< qM^M I e-'"'/" I QM-i^M-i > ■ • ■< qiai \ e-'"''^ \ qoPo > (8) 

where 

qm = X, = a*, go = x', a© = /3 . (9) 
Now insert M completeness relations in the momentum representation 



j dp \ p >< p 1= 1 



-5{/-^}n{/-} 



in (8) in order to obtain 

k=i y." } k=i 

< qnocn I e-^-^'/^ I PnoIn-i >< Pn \ Qn-i > 

■ ■ ■< qidi I e-'"'l^ I Pido >< Pi I go > (10) 

Now we will take the limit that M — > oo, e ^ but with t — (M — l)e being finite. 
For small e we should expand the exponential in (10) to first order in e and write 

< Qkdk I e"'^^/^ I Pkdk-i >~< Qk I Pk >< ttfe I ttfe-i > 

exp^-^eH{qk,Pk,dl,dk-i)^ (11) 

where 
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< a I afc-i > 
Using the overlaping relations 

< a \ f3 >= exp {a* ■ (3 — 



<P\<l>=^exp[-^Pg] 



one obtains 



27Tn \J 2'Kh TT^ 

fe=l 



fc=l ^ 



2 

Now one has to integrate over Pfc. So, using (1) we must evaluate 



dPk \ ie 



2M 



where 



o 



^ / ' ]\/r dnmC^n^rn 

n,m=l ^^^O 

This allows one to rewrite (14) in the standart form [|l 



Mo \''^ [iMoe (.{qk-qk-i] 
exp < /I + 



,2niheJ \ 2h \ e J 

Now, substituting (16) in (13) and taking the limit of M ^ oo and 



|3| 1/51 

= / Dq l D'^ae 2 2" exp{S'[g, a]} 

Jx' J /3 
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where 



k 

k=l 



M-1 / •• 

D^a= lim TT 



k=l \ 



with 

S = I dt' I (g + (5))^ + - fcJo; — a* ex) — —Hr^q* 

JO \ fhZi Zi ft 

Finally, using (3) and (15) one reaches 

Now, substituting (21) in (17) and (6) we get the result (4.12). 
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